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Abstract 

We present a quenched lattice calculation of the nucleon isovector vector and axial- vector charges 
g v and g A . The chiral symmetry of domain wall fermions makes the calculation of the nucleon 
axial charge particularly easy since the Ward-Takahashi identity requires the vector and axial- 
\Q ' vector currents to have the same renormalization, up to lattice spacing errors of order 0(a 2 ). The 

o ' 

DBW2 gauge action provides enhancement of the good chiral symmetry properties of domain wall 



fermions at larger lattice spacing than the conventional Wilson gauge action. Taking advantage 
of these methods and performing a high statistics simulation, we find a significant finite volume 
effect between the nucleon axial charges calculated on lattices with (1.2 fm) 3 and (2.4 fm) 3 volumes 



(a rs 0.15 fm). On the large volume we find g A = 1.212 ± 0.027(stat) ± 0.024(norm). The quoted 
systematic error is the dominant (known) one, corresponding to current renormalization. We 
discuss other possible remaining sources of error. This theoretical first principles calculation, 
which does not yet include isospin breaking effects, yields a value of g A only a little bit below the 
experimental one, 1.2670 ± 0.0030. 

PACS numbers: 11.15.Ha, 11.30.Rd, 12.38.Aw, 12.38.-t 12.38.Gc 
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I. INTRODUCTION 



The axial charge g A of the nucleoli, or more precisely its ratio to the vector charge, g v , 



appears to be a good test of our understanding of the structure of the nucleon. 



r irst of 

H 



all, it is very accurately measured from neutron (3 decay, g A /g v = 1.2670 ± 0.0030 
And, among the nucleon form factors or moments of structure functions, it is technically 
the simplest from the point of view of a lattice QCD numerical calculation. 

Four form factors appear in neutron (3 decay: the vector and induced tensor form factors 
from the vector current, 

{p\V+(x)\n) = u^^q 2 ) - q X (Tx^g T (q 2 )}u n e- iq - x , (1) 

and the axial- vector and induced pseudo-scalar form factors from the axial- vector current, 

(j)\A+(x)\n) = u P [l^ 5 g A (q 2 ) - iq^ 5 g P (q 2 )}u n e- iq - x . (2) 

Here q = p n — p p is the momentum transfer between the proton (p) and neutron (n). In the 
limit \q\ — > 0, the momentum transfer should be small because the mass difference of the 
neutron and proton is only about 1.3 MeV. This makes the limit q 2 —>■ 0, where the vector 
and axial- vector form factors dominate, a good approximation. Their values in this limit are 
called the vector and axial charges of the nucleon: g v = g v (q 2 = 0) and g A = g A (q 2 = 0). 
Experimentally, g v = cos9 c (with the Cabibbo mixing angle 9c), and g A = 1.2670(30) x g v . 

Since they are defined at zero momentum transfer, a naive expectation is that g v and g A 
are easier to calculate on the lattice than form factors which require non-zero momentum 
transfer. Despite this, quenched QCD lattice calculations with Wilson fermions at finite 
lattice cutoff (a -1 ~ 2 GeV) have underestimated g A by about 20% 0,0,0] (see Table H]for 
a summary of previous calculations). This suggests systematic errors, which may arise from 
(1) the quenched approximation, (2) operator renormalization, (3) non-zero-lattice-spacing a 
and loss of chiral symmetry for Wilson and Kogut-Susskind fermions, and (4) finite volume, 
remain in the lattice calculation. 

The first three errors have been addressed in previous calculations. The SESAM and 
LHPC collaborations found that unquenching does not solve the problem as the estimated 
value g A decreases by 5-10% 0, Q]. On the other hand, reducing the lattice spacing error 
seems to increase the value, but only by a small amount, ^ 5% 0, 0]. Perhaps more 
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important is the calculation of the renormalization factor Z A for the axial current. The one- 
loop perturbative renormalization factor, used in the case of Wilson fermions [3, S, U, 0, 0| , 
was probably overestimated. The QCDSF-UKQCD collaboration reported that the non- 
perturbatively calculated renormalization factor (2'° on P ert ~ 0.8) is roughly 10% smaller than 
the one-loop one [Z^ ert ~ 0.9) in the case of the non-perturbatively 0(a) improved Wilson 
fermions [8] at a -1 ~ 2 — 3 GeV. Thus, the systematic error in the determination of the 
renormalization factor appears to be more important than the first two effects mentioned. 
The first two systematic errors listed above likely cannot resolve the issue that previous 
lattice calculations of g A underestimate the experimental value. 

The loss of chiral symmetry on the lattice is potentially significant. As is well known, 
g A /g v = 1 in the absence of chiral symmetry breaking in QCD. Further, in the realistic case 
of spontaneously broken chiral symmetry, the ratio is still constrained by the axial Ward- 
Takahashi identity; d^A^x) = 2mP a (x). The Goldberger-Treiman relation derives from 
the nucleon matrix elements of the currents on both sides of this identity in the soft pion 
limit 9}. We can easily understand the deviation of the ratio from unity in the context of 
the Gell-Mann-Oakes-Renner relation [10] which is also related to the axial Ward-Takahashi 
identity. Thus, the explicit breaking of chiral symmetry at non-zero lattice spacing a for 
Wilson fermions may induce significant errors which are only removed in the continuum 
limit. 

In this work we use domain wall fermions fDWF), a fermion discretization scheme with 
almost perfectly preserved chiral symmetry JjJ, [12J, [l3j . This scheme introduces a fictitious 



fifth dimension in addition to the four dimensions of space-time. In the limit where the fi 
dimensional extent L s is taken to oo, DWF preserve the axial Ward Takahashi identity 



Ih- 
ljl 



at non-zero lattice spacing. With finite L s the suppression of explicit chiral symmetry 
breakin g is effectively exponential in quenched simulations if the gauge field is sufficiently 

19, 20]. This is always true if the lattice spacing is sufficiently small. 



smooth 15lllfi. 



17. 



In low energy cases like the one investigated here, the small breaking of the symmetry at 
finite L s is parametrized by a single universal "residual mass" parameter, m res , acting as an 



18, 



2l|. 



additive quark mass and which is defined from the axial Ward-Takahashi identity 
Furthermore, the DWF scheme greatly simplifies the non-perturbative determination of the 
renormalization of quark bilinear currents j^j. For example the renormalization factor of 



local vector and axial- vector current operators should be equal, Z A = Z v |22j. This means 
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the ratio of the nucleoli axial and vector char ges calculated on the lattice directly yields the 
continuum value, i.e. it is not renormalized J23I- 1^. By employing the DWF scheme, the 
ambiguity in the renormalization of quark currents which may be present and problematic in 
other fermion discretization schemes is eliminated. We emphasize that the DWF calculation 



of the nucleon axial charge should not suffer 
renormalization and loss of chiral symmetry 



rom the systematic errors due to the operator 
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24 1 . 



However, as is described in more detail in section III, in our first DWF calculation with 
the single-plaquette Wilson gauge action at (3 = 6.0 and lattice volume 16 3 x 32 x 16 (which 
correspond to a -1 ~ 2 GeV and spacial volume ~ (1.6 fm) 3 ), we found that g A exhibits a 



fairly strong dependence on the quark mass 



24j . A simple linear extrapolation of g A to the 



chiral limit yielded a value that was almost a factor of two smaller than the experiment 



241- 



This implied the presence of a large finite volume effect. To our surprise, we found no 



systematic study of such an effect in t 
dependence in the naive quark model [2 



:he literature. Note also that there is no volume 
nor in the MIT bag model j^l . In the former the 
ratio is determined by a simple spin-isospin algebra, and in the latter it arises from a simple 
overlap integral of the upper and lower component of the bag Dirac wave function. 

To address the finite volume issue we need to have at the same time a sufficiently high 
lattice cutoff to preserve chiral symmetry reasonably well and at least two lattice volumes, 
preferably ones that are large compared to the charge radius of the proton. The Wilson gauge 
action will not work for this purpose since the chiral symmetry of DWF in the quenched case 
degrades rapidly as lattice spacing a increases, while the computational cost necessitated 
by a very large lattice volume would be prohibitive. Fortunately various "renormalization- 
group-inspired" improved gauge actions preserve the chiral symmetry of DWF well while 
not demanding a large cutoff jigL l^. Thus both requirements, chiral symmetry and large 
physical volume, can be met at reasonable computational cost. Of the relatively well- 
established candidates in this class of improved gauge actions, we choose the "doubly-blocked 



Wilson 2 (DBW2)" action |20j,|22|. 



The rest of this paper is organized as follows: in section II the lattice method for cal- 
culating g A /g v is described. In section III the numerical results obtained for both Wilson 
and DBW2 actions are described in detail. Finally, in section IV we summarize the present 
work and discuss future directions. 
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II. GENERAL ANALYTIC FRAMEWORK 



A. The vector and axial charges 

As mentioned in the introduction, four form factors are needed to describe neutron (3 
decay: the vector and induced tensor form factors for the vector current, 

(p\ V ^( Q )\ n ) = Up^^gyiq 2 ) - q\(rx^g T (q 2 )]u nj (3) 

and the axial and induced pseudo-scalar for the axial current, 

(p|Aj(0)|n> = u p [-f^ 5 g A (q 2 ) - iq^g P {q 2 )]u n . (4) 

The right hand side of each is the most general form consistent with Lorentz covariance. 
The momentum transfer q = p n — p p becomes very small in the forward limit because of 
the small mass difference between the neutron and proton. In the limit q 2 — > 0, which we 
take in this work, the vector and axial form factors dominate the matrix elements. We 
are neglecting the mass difference of the neutron and proton, and hence that of up and 
down quarks (we also neglect the electromagnetic mass difference.) For zero quark mass 
m = m u = rrid = the action is symmetric under global chiral 577(2) x £77(2) flavor 
rotations acting on the quark fields. If m ^ 0, the symmetry is broken down to the vector 
(flavor) 577(2) sub-group, and the associated vector charge, g v is still conserved (g v = 1). 
This situation is sometimes called CVC, conserved vector current. In the real world even 
this symmetry is softly broken by the small mass difference between up and down quarks, 
m u — irid. The explicit violation of the axial- vector symmetry by non-zero quark mass is 
sometimes called PCAC, or partially conserved axial-vector current. As is well known the 
axial SU (2) symmetry is also spontaneously broken. Thus, the axial charge may in general 
deviate from unity, g& ^ 1. 

If the vector symmetry is preserved, a simple exercise in SU(2) Lie algebra leads to the 
following (see Appendix lA"|) : 

(p\A+\n) = 2{p|4b) = (p\A;\p) - (plAjlp) (5) 

where = uj^d, A 3 ^ = \{At — Afy, A^ = uj^u and A^ = d^^^d. u and d stand for 
the up and down quark fields. A similar relation holds for the vector case, 

(p\V;\n) = 2(p\V'\p) = (p\V;\p) - {p\V*\p) (6) 
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where = u^^d, V£ = \(V™ — V*), V™ = uj^u and Vjf = d^y^d. The isovector vector 
charge g v and the isovector axial charge g A are defined by the strength of the right-hand side 
of Eqs. EJandiniin the forward limit (q 2 — > 0). In addition, the polarized quark distributions 
in the proton for each flavor /, Aipf, are defined by the forward matrix element of the flavor 
axial- vector currents Af,: 

(k, s 1^/7^75^/1 k, s) = Aip f u p (k, s)-/^ 5 u p (k, s) (7) 

where k and s are proton four momentum and polarization. From CVC we find the relation 
g A = Au — Ad. 

Now consider the conserved electromagnetic current expressed in terms of the flavor 
vector currents V.{ : 

3r=Y,Qf v i=i v ?-3 v i + -~ ■ ^ 

Here Qf denotes the charge (in units of proton charge e) for a quark of flavor /, and the 
ellipsis denote possible flavors of heavier quarks which we henceforth ignore. Since the 
corresponding electromagnetic U(l) gauge symmetry assures conservation of electric charge, 
for the neutron we find 

lim <n|#» = 0. (9) 

It follows that 

lim (n\V*\n) = 2 x lim (n\V^\n). (10) 
On the other hand, under the assumption of CVC we have the following: 

(n|V» = (pK d \ P ) (11) 
(n\V«\n) = (p\V;\p). (12) 

Thus we reach the following relation: 

lim (p\j™\p) = lim (p|V» = lim <p\V? - V*]p). (13) 

q 2 ^0 q 2 ^0 q 1 — >0 

Likewise, it follows that the vector charge, g v , must be unity (in units of cos 6*c and e) under 
CVC since the proton electric charge is unity. As already mentioned, we expect a very small 
breaking from CVC because of the physical up and down quark mass difference. In the axial 
case, we expect non-conservation of g A due to the small but non-zero up and down quark 
masses as well as the spontaneous breakdown of chiral symmetry. 



B. Nucleon matrix elements 



In this subsection we describe our method of lattice numerical calculation of the axial 
and vector charges of the nucleon. Hadronic matrix elements calculated on the lattice are 
determined from ratios of the relevant three-point to two-point correlation functions. Since 
the charges are defined at zero-momentum transfer, we do not have to introduce non-zero 
momentum projection for the nucleon source and sink for these correlation functions, nor for 
the current insertion. On the other hand since we are dealing with a spin-1/2 baryon, both 
correlation functions possess non-trivial Dirac spinor structure, so appropriate projections 
are necessary. 

The zero-momentum two-point function for the nucleon is given by the sum over all 
spatial coordinates, x, 

(Af(tW(0)) aP = ^(0\K&tWp@,0)\0) (14) 

X 

where J\f(x, t) can be any operator with the same quantum numbers as the nucleon, namely 
unit baryon number, J p = (l/2) + , and isospin doublet, a and (3 denote Dirac indices. 
Color and flavor indices are suppressed in the following unless noted otherwise. The two- 
point correlation has the asymptotic form 

(Af(tW(0))*P = 4f [1 + sgn(t) l4 Ue- M ^ (15) 

at large Euclidean time, t. Here Mjq denotes the ground state mass of the nucleon. The am- 
plitude A N is defined as (0|A/"(0) |iV) = A N u^- In general, the baryon two-point function 
receives contribution from both positive and negative-parity states. By taking the trace with 
a projection operator P + — (1 + 74)/2, we eliminate contributions from the opposite-parity 
state in forward time direction. Details of the parity projection are described in Let 
us abbreviate the notation for the two-point function of the particle contribution from the 
desired (positive-parity) state as 



G N (t)= 1 -Tr[P + (Af(tWm]- (16) 

The factor of 1/4 is our choice of normalization. At large t this asymptotically approaches 
a simple exponential, 

G N (t)~exp(-M N t). (17) 



For the proton, a standard choice for the interpolating operator is 

Af(x) = e ahc [u T a (x)C^d b (x)]u c (x) 



(18) 



where C is the charge conjugation matrix defined as C = 7472, abc the color indices, u and 
d the up and down quark fields. 

Next, let us define the zero-momentum three-point correlation function for quark bilin- 

ears, { T f \x) = ^ f {x)Y-ij) f {x)\ 

m)Or(tW(o)u = EE<°i T {^ *w»(o, °)}i°> ( 19 ) 

x x' 

where T is any of the sixteen possible matrices in the Clifford algebra defined by the Dirac 
gamma matrices. When t ^> t' ^> 0, the particle contribution of the zero- momentum three- 
point function becomes 

(M(t)O r (t'W(0)) a p - A N lim g r (q 2 ) exp(-M N t)(P + TP + )^. (20) 

Note two important points: first, the three-point function vanishes for T other than 1, 74, 
7«75 (i = 1, 2, 3), and 07,- (i, j = 1, 2, 3) because P+rP + = for Fs that do not commute 
with 74. Second, the r.h.s. of the above asymptotic formula does not depend on the insertion 
point of the operator Or- Any t'-dependence arises from excited state contamination, i.e. 
away from the asymptotic regime. 

In this paper, we calculate the isovector (quark-flavored) vector charge g v and the isovec- 
tor axial charge g A of the nucleon. We define the spin projected three-point function 
for the relevant components of the vector current V* = V'/7/*V'/ an< ^ the ax i a l current 
Aft = ^Pflnlb^f by taking traces with the projection operators Pp = P+r^ 1 : 

G f v (t,t') = ^Tr[P v (Af(t)V 4 f (t')IT(0))} (21) 

G f A (t,f) = InPAAmMit'wm] m) 

where P v = P+ and P Ai = P+7i75 (i = 1, 2, 3). In order to extract g T (Y is either V or A) on 
the lattice, we have to identify a plateau in the ratio of the three- and two-point functions, 



lattice G^t')-G%t') 



G N (t) 

in the range of t > t' with fixed t = i s i n k — Source- 



(23) 



In general, lattice operators receive finite renormalizations relative to their continuum 
counterparts since the exact symmetries of the continuum are usually realized only in the 
continuum limit, a — > 0. Thus 

g 1 ™ = Z r g l ; ttice (24) 

requires some independent estimation of Z r , the renormalization of the quark bilinear cur- 
rents, 

[^r^] rcn = z r [#v>] lattice . (25) 

Since DWF possess full chiral symmetry at non-zero lattice spacing, a lattice conserved 
vector current and partially-conserved axial-vector current A^ which receive no lattice 
renormalization can be defined, namely Z V =Z A =1 However these conserved currents 
are point-split and require sums over the extra 5th dimension of DWF, so they are somewhat 
costly to work with in practice. Alternatively the local currents and which are 
naive transcriptions of the continuum operators are easier to deal with but receive a finite 
renormalization since they do not correspond to an exact symmetry of the action. However, 
the Ward-Takahashi identity satisfied by both types of currents is enough to ensure that the 
lattice renormalizations of the local currents are equal, Z v = Z A , up to terms of order 0(a 2 ) 
in the chiral limit and neglecting explicit chiral symmetry breaking for DWF at finite L s . 

Note that the vector charge computed from the local current provides an independent es- 
timate of Z v since the renormalization of an operator does not depend on any particular 
matrix element and the renormalized, or physical, value of g v is 1 by CVC. 

7 (07) 

v G» (t,t>) - G d (t,t'Y { } 

V V 

Comparison of Z v thus obtained to the value of Z A from the relation [is| . 

(A(*75^](0)> = Z A (A»(t)[^}(0)), (28) 

yields an estimate of the 0(a 2 ) systematic errors arising from the method described here. 
These are discussed in section III. 

Next we describe the particular interpolating operators, or quark sources and sinks, used 
to calculate lattice correlation functions. In our earlier work we used so-called wall-wall 
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correlation functions constructed with quark sources generated from a unit source at each 
spatial site on a fixed source time-slice and summed over all spatial sites at the sink time- 
slice. Since the wall source or sink is gauge variant, we fix to the Coulmob gauge. Later we 
switched to wall-point correlation functions since they yield smaller statistical errors. For 
three point functions this approach is implemented with a sequential source. We discuss 
both types of correlation functions in turn. 

First, we introduce the forward quark propagators from a wall source to a wall sink and 
from a wall source to a point sink, which may be written with the gauge fixed point-to-point 
quark propagator S pp : 

S ww (t,t')% = Y^S pw {x^t')%, (29) 

X 

S p Jx, t; t')% = S PP & *5 1, 0$ (30) 
y 

where the subscripts and superscripts denote Dirac and color indices, respectively. The 
quark three-point function resulting from insertion of the quark bilinear operator ipfTipf is 
defined as 



x,y,z 



= E -rwSlw & t)^Mr) lY s p jx, f ; t'X>p (si) 

z 

where the second line results from S (x,t;y,t') = j^Sp (y, f; x, t)j 5 . Thus, S r (t,t',t") 
is constructed by combining wall-to-point quark propagators generated from two different 
source time-slices t and t" at either end of the lattice with the operator inserted in between 
them. 

The two-point function for the nucleon in Eq. (jl4j) is expressed in terms of quark propa- 
gators as 

(AT(tW(0)) aa , = e^ a , vd {Ci h )frS^(t,U)%{Ci h )w (32) 

Following Ref. the three-point function in Eq. (jl9j) is easily obtained from the two-point 
function by replacing the ordinary quark propagator by the operator inserted one, <S r (t, t', 0). 
Inserting the d and u quark currents, we obtain 

(Af(t)o¥\t'W(0)) aa , = e abc e a , b ,ACj 5 )^ d \t,t>,0)^,(C l5 )p, Y (33) 
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X 



and 



(AT(t)4 u) (t')AT(0)) ( 



abc ea'vAC l5 )^S%l(t, 0)%,(Cy 5 )w (34) 

f , o)^sJt ( t , oft, + f , o)^£t(t, oft 

+<S£t (*> 0)^^(t, f , Oft, + S%l(t, 0)^5 f , oft/ . 

The nucleon three-point function is the sum of the up and down quark contributions. The 
spin projected three-point functions are obtained from EqsJ21land 1221 

To enhance the signal, a point sink is more desirable than an extended sink. The wall- 
point type of three-point functions is implemented using the so-called sequential source 
method 0, 0, 0. In addition we use box source instead of wall source to enhance the 
coupling to the ground state: 



o<y<B 



(35) 



We adjust the box size B to about 1 fm. In describing the construction of the sequential 
source, it is convenient to introduce the "diquark" propagators: 



T> aa ,(y,t;O)% a = s abc e alVc/ (Cj 5 ) d p(Cj, 



'5)S'/3' 



S { : ] M *5 0)fpS<£(y, t; 0)# + S^(y, t; 0)%,S^(y, t; oft] (36) 



and 



U aa ,(y,t;0)fy = e^eawS^iy^O)^, (37) 
x \c l5 )^(C l5 )pyS^(y, t; 0)™', + (C l5 ) 5 ~ ( (C l5 ) s ,yS^(y, t; 0)^6^6^ 
+ (C^)p 7 (C l5 ) 5 ^S^(y, t; 0)« a X>P> + (Cl 5 )s,(C l5 )^ Y S^(y, t; 0)^<W 

The "down diquark" (T>) and "up diquark" (U) are defined by the down quark removed 
propagator and one up quark removed propagator from the nucleon two-point function. 
Now using the diquark we can reconstruct the point-to-wall type of the nucleon two-point 
function as 

y y 

(38) 
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In terms of the diquarks the three-point functions of an arbitrary quark bilinear operator 
i^/Ttpf at a location (z,t') can be written for the down quark 

(AT(t)4 d) (t')AT(0))^ = ^X>^(y,t;0)^^(y,t;z,Or*(n«'^(^^°)^ ( 39 ) 
and for the up quark 

{N(t)OP{t')Jf(V)) aa , = ^ ^(^, t; 0)J^W(y, t; t')^(r) w ^(^ f'; 0)^,. (40) 

y,z 

For the construction of the three point functions we need the backward propagators from 
the sink point (y, t) to the operator insertion point (z, t') . However, it is highly expen- 
sive to prepare the required point-to-point quark propagators from all points (y, t) . This 
difficulty is easily circumvented by directly computing the generalized quark propagators 
D aa » (y, t; 0) (y, t; z, t') and tl aa < (y, t; 0) S 1 ^ (y, t; z, t') with the sequential source method. 

Before describing details of the sequential source propagator, we should apply the spin 
projection Pp to diquarks in order to reduce the cost from having to calculate all 4 x 4 
matrices for external spinor indices (a, a'). In this article, we only need two kinds of spin 
projections, i.e. Py and Pa 3 so that it reduces the amount of calculations by a factor of 
eight in comparison with the unprojected case. The spin projected source for a down quark 
insertion is 



*d(i7,t;0)$ = -2V(y,*;0ft(P r ) yT , (41) 



1, 

4' 

and for the up quark 

*«(£*; o)$ = jUrf(v,t;0)%mfv ( 42 ) 

Finally, the sequential source down quark propagator is 

V d (t,0;z,t'yi= [X;75^5M(£^;^^75,y^**S(i/»^0)S] » ( 43 ) 
\y,t" J 

and the sequential source up quark propagator is 

E tt (f , 0; z, t)% = I 75^M(£, t ' ; £ f)^,^**^, t"; 0)5 I • (44) 
\y,t" / 

which may be calculated by solving the matrix equations 

E S /(*' f » OS Ai f (f , t"; z, t')% = 5 tt ^f(z, t>; 0)£ (45) 

x,t" 
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where is the Dirac matrix. Consequently, in terms of the sequential source propagator, 
the spin projected three point function for the down quark is written 

G d r (t,t') =^S^ ) 0;£,Oa(r)/*r5{2(^* / ;0ft (46) 

2 

and for the up quark is 

G^(M') = E^^^^SCrK^^^Oft. (47) 

2 

In the case of keeping the up and down quark masses equal, the total cost for computing 
the sequential source propagator is a factor of two over the cost for wall-wall correlation 
functions. However, the resulting box-point correlation functions yield smaller statistical 
errors. 



III. NUMERICAL RESULTS 



We have performed quenched lattice calculations using two different gauge actions, the 
standard Wilson and the improved DBW2 j^. Details and some relevant results of both 
simulations are summarized in Tables II, and IIV1 We describe the nucleon matrix 
element results for each one separately, then compare them and draw some conclusions. 



A. Wilson gauge action results at (3 = 6.0 

We have performed a quenched simulation on a 16 3 x 32 lattice with the standard single- 
plaquette Wilson action at /3 = Q/g 2 = 6.0 which corresponds to a lattice cut-off of a -1 = 
1.922 GeV set by the p mass jl8j ]. Quark propagators were generated with four bare masses, 
m f =0.02, 0.03, 0.04 and 0.05, using DWF with L s = 16 and M 5 = 1.8. The nucleon matrix 
elements were averaged on a set of 400 gauge configurations. Hadron masses computed on 
these lattices are tabulated in Table Preliminary results for the nucleon charges were 



first reported in 24 1. 



We calculated wall-source quark propagators on each Coulomb-gauge-fixed configuration 
for both periodic and anti-periodic boundary conditions in the time direction for the quarks. 
A simple linear combination of these propagators then yields a forward (or backward) in 
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time propagator. To compute the correlation functions, we employed the wall-wall method 
described in the previous section with source locations fixed at t src = 5 and t' src = 21. 

In Figure H we show the dependence of the vector renormalization, Z v = 1/ g^ ttlce on the 
location of the current insertion. A good plateau is observed in the middle region between 
the source and sink. The quoted errors are estimated by a single elimination jack-knife 
method. The dashed lines represent the average value and statistical error in the time-slice 
range 5 < t — t src < 11. The mass dependence of Z v is rather mild as seen in Figure 121 and 
given in Table IVT1 The values 0.7601(31) for a linear fit and 0.7610 (52) for a quadratic fit 
at m/ = agree well with Z A = 0.7555(3) Q], which was obtained from a calculation of 
meson two-point correlation functions. The discrepancy AZ = |1 — Z A /Z v \ is less than 0.6% 
which implies the 0(a 2 ) error that remains after taking the rrif —>■ limit is quite small. 

As is seen in Figure El plateaus are evident for the spin-dependent distribution functions, 
Am and Ad, in the range 5 < t — t STC < 11. Thus, we compute the charge ratios {g A /g v ) la,ttlce 
at each my by taking a weighted average over this time slice range. In Figure 0] a strong 
dependence on rrif appears. A simple linear extrapolation to rrif = yields 0.812(112), 
which is roughly 2/3 of the experimental value. However, a simple linear ansatz may not 
describe the data which show increasing downward curvature for lighter quark mass (note 
that the points are correlated in this quenched calculation since they are computed on the 
same gauge configurations). In general chiral logarithms may appear and were considered. 
In fact, the data are not compelling for such terms, arising in either quenched or full chiral 
perturbation theory The results for each mass are reproduced in Table IVI1 

This implies the existence of other systematic errors. As was mentioned in the introduc- 
tion, a large systematic error in previous lattice calculations of g A came from the determi- 
nation of the renormalization constant Z A . As shown above using DWF, the value of g A is 
determined in a fully nonperturbative way, with or without explicit renormalization. The 
systematic error stemming from the incomplete cancellation of renormalization factors in the 
ratio is less than 1% as we saw by comparing Z v = 1/ g^ ttlce and Z A calculated from meson 
two-point functions. In addition, comparing the chirally extrapolated values of {g A / g v ) Uttlcc 
and Z A g l ^ ttlce leads to an even smaller error, although it relies on the linear extrapolation 
which was not very compelling. Another possible systematic error is the contribution of 
excited states, the presence or absence of which was checked by slightly enlarging the sep- 
aration between wall sources, t S rc(= 5) and t S i n k(= 27). While the larger separation induces 
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more noise in the signal, the central value of g A is essentially unchanged for each quark mass; 
thus we cannot detect a systematic effect outside of the statistical errors. Still, this source 
of error appears to be small. 

Detailed detection of quenching effects: quenched chiral logarithms, unsuppressed 
fermionic zero modes, and the absence of the physical pion cloud, is beyond our scope 
at present since these require very light quark masses and correspondingly large statistics. 
Thus, by process of elimination we are lead to focus on finite volume effects which we discuss 
in the next section. The volume employed for the calculations in this subsection is roughly 
(1.5 — 1.6 fm) 3 which can barely accommodate a proton with mean square radius estimated 
to be about 0.8 fm ill. 



B. DBW2 action results at (3 = 0.87 

To determine g A in a large physical volume, say <^ (2fm) 3 , we have performed a DWF 
simulation on a lattice with larger spacing. In general, it is difficult to maintain the good 



chira 
tion 



properties of DWF as a increases at fixed L s , especially with the Wilson gauge ac- 



18, 



191 ] . It has been shown that the Iwasaki gau ge a ction enables studies of quenched 



19L |35[ . Recent quenched studies by 



DWF with smaller L s than the Wilson gauge action 
the RBC collaboration have shown that the chiral symmetry of DWF are even better with 
a similar type of renormalization group improved gauge action, DBW2 j^J. The chiral 
symmetry of DWF with DBW2 is significantly improved over the Iwasaki action. A very 
small additive quark mass m res ~ 0.8 MeV is achieved on a lattice with a -1 « 1.3 GeV and 
L s = 16. Good scaling behavior of the light hadron spectrum is observed as well j^J. 

To study finite volume effects numerical simulations were performed at (3 — 0.87 (a ~ 0.15 
fm) on two lattice sizes, 8 3 x 24 and 16 3 x 32, with L s = 16 and M 5 = 1.8. Our results 
are analyzed on 400 quenched gauge configurations for the smaller lattice (La ~ 1.2 fm) 
and 416 configurations for the larger lattice (La ~ 2.4 fm). Hadron masses computed in 
this calculation are summarized in Table IVII1 Meson masses (m, and m p ) for the 16 3 x 32 
lattice are evaluated from 100 configurations. 

In this calculation, we utilize the sequential quark propagator method to compute three- 
point functions as described in section II. We checked for consistency with the wall-wall 
method on the smaller 8 3 x 24 lattice. The sequential- and wall-type quark propagators in 
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Coulomb gauge were computed at five evenly spaced values of rrif ranging from 0.02 to 0.10. 
The smallest quark mass corresponds to a pion mass m n ~ 390 MeV. The Nucleon source 
and sink were separated by about 1.5 fm, which corresponds to the same physical separation 
in time used in the calculation with the Wilson gauge action at /3 = 6.0. A preliminary 
version of the results presented below was first reported in j^, [mJ . 

First, we check whether Z v = Z A is true even on this coarse lattice. The vector renormal- 
ization 1 / ^f^ ttlCG is plotted against the location of current insertions in Figure The data 
are calculated on the larger spatial volume with the sequential quark propagator method. 
We take a weighted average of i/g l ^ ttlce with the three middle points (t — t STC = 4,5,6) to 
evaluate the vector renormalization Z v . The dependence of Z v on mj is shown in Figure |H1 
and given in Table ITX1 In general = Z v V fl + 0(a 2 ,mfa 2 ) and A^ = Z A A^ + 0(a 2 ,rrifa 2 ), 
where V M and A^ denote the conserved vector currents. It is not so apparent in our data. 
A linear extrapolation yields Z v = 0.7952(13) at rrif = 0, while a linear plus quadratic 
extrapolation gives the value 0.7991(25). The RBC collaboration obtained the renormal- 
ization factor of the axial- vector current Z , nonperturbatively from a calculation of meson 

nn 

two-point correlation functions [la, |2(J] . It was found Z A = 0.77759(45) in the massless 
limit (stj which is smaller than the value of Z v obtained above by 2-3%. This discrepancy 
may be caused by an order 0(a 2 ) lattice artifact. 

To explore this possibility further, we evaluate the renormalization factor of different 
vector currents. According to Secjnl the conserved current V M guarantees that nucleon 
matrix elements of j^ 71 , and — should be identical. The local lattice currents 
are renormalized as = Z^Vj + 0(a 2 ) = Z v Vj[ + 0{a 2 ) in the chiral limit. FigureEl 
shows the values of Z^ as well as the value of Z A from 20]. The difference among values 
of Z v appears independent of rrif within statistical errors, and the discrepancy between the 
smallest and the largest is comparable to that between Z v and Z A noted above. We also 
note that the discrepancy is larger at this lattice spacing, by roughly a factor (1. 922/1. 3) 2 , 
than the corresponding one for the Wilson gauge action results discussed earlier. Of course, 
since the two gauge actions have different 0(a 2 ) errors, the comparison is only a crude one. 
We conclude that Z v = Z A is satisfied up to small discretization errors of 0(a 2 ) on this 
coarse lattice. 

In Figure [HI we plot unrenormalized spin-dependent densities Am and Ad, which are 
calculated with the sequential source propagator j^], against the location of the current 



16 



insertion. In this calculation, the sequential source propagator was calculated with a box 
source and a point sink, so the resulting three-point function has no time reflection symmetry 
about the midpoint between t STC and t sin ^ because excited state contamination is worse for 
the nucleon propagating between the operator and the point sink. In Figure |H1 the plateaus 
appear shifted toward the wall source, as expected. Next we evaluate the bare value of 
^lattice _ ^ u _ a £ gg^h m ^ shown in Figure El for both lattice volumes and tabulated 
in Tables IVHHUXI g^ ttlce evaluated on the smaller volume is clearly smaller for each value 
of m/, and the difference increases as mj decreases. In contrast, g v does not show much 
dependence on the volume. 

To compare to our previous DWF results with the Wilson gauge action, we plot the 
value of {g A / g v ) lattlce as a function of (m n /m p ) 2 in Figure ITU1 The smaller volume results 
using the DBW2 gauge action {(3 = 0.87) are the same (within statistical errors) as our 
previous results using the Wilson gauge action (/? = 6.0) on a slightly larger volume. The 
large volume DBW2 results exhibit mild quark mass dependence while both smaller volume 
results show a marked decrease toward the chiral limit. We conclude that our previous 
DWF- Wilson-gauge-action results were significantly adversely affected by finite volume. 

Finally, we extrapolate g A en to the chiral limit. For this purpose, we have two methods. 
One is to extrapolate the charge ratios {g A / g v ) lattlce to the chiral limit where the relation 
Z v = Z „ is valid. The second method is the conventional one utilized in all other calcula- 



tions [2j, y, |4], |5j, |6j, 13, |8j . The chiral extrapolation is performed on g lattlce x Z A . Recall that 



A 



the latter requires the value of Z A , whether nonperturbatively or perturbatively calculated, 
while the former does not. In the present case, we use the nonperturbative value of Z A 
from Q. 

We plot {g A / g v ) lattlce and Z A x g^ ttlcc together in Figure ^2 and perform a simple linear 
extrapolation in each case. The two methods provide consistent results in the chiral limit: 
the ratio method gives g™ n = 1.212(27) while the conventional method gives g™ n = 1.188(25). 
In light of our earlier discussion, the systematic difference if there is one, is related to our 
choice of renormalization. A two percent error stemming from Z v ^ Z A yields 0.024. This 
is also the difference in the central values just obtained. Thus, we quote 

g ren = 12U ± .027(stat) ± 0.024(norm) (48) 
which underestimates the experimental value of 1.267 by less than five percent. We have 



17 



not attempted to estimate residual non-zero lattice spacing, finite volume, explicit chiral 
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symmetry breaking, and quenching effects. The first three are probably small 
The only remaining error not under good control is the quenching one which does not appear 
to be large 0, 0], in light of the relatively good agreement with experiment shown above. 
This view does not change unless significant non-analytic behavior, which we did not detect 
here, arises near the chiral limit. 

We note Jaffe recently showed that in the chiral limit the nucleon axial charge is delo- 
calized, and he argued this leads to a large reduction in g A calculated in a finite volume 
surrounding the nucleon Subsequently, Cohen showed that in a finite volume with 

periodic boundary conditions pertaining to lattice calculations, this phenomenon does not 
lead to a reduction in g^ j^II- However, as emphasized in 37 1, this does not preclude other 
large finite volume effects. 

As mentioned above, this calculation of g A is performed for relatively heavy quark masses; 
the quenching error at this unphysically large mass scale is probably small. However, one 
may worry that such a calculation does not capture relevant physics in the region where the 
quark mass is much lighter, and the so-called "pion cloud" surrounding the nucleon becomes 
important. Nevertheless the values of g A en at these heavier quark masses already lie just a 
few percent below the experimental value and show little dependence on the quark mass. 
This presents an important question concerning the role of the pion cloud: is it a few percent 
effect, as seems plausible from our first principles calculation, or is it larger, as estimated 
from phenomenological models 

The dependence of the product m N g A on the lattice volume is of interest (See Fig. IT2|) . 
While the smaller volume results always lie below the larger volume ones, within one stan- 
dard deviation they almost always agree. There is only one exception at arrif = 0.08 in 
the bare lattice result. No volume dependence is detected. This is in clear contrast to 
the situation of the axial charge alone. Since the product is the one that appears in the 
Goldberger-Treiman relation, it would be interesting to see how its counterpart, the induced 
pseudo-scalar form factor, behaves at small momentum transfer. 
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IV. CONCLUSIONS 



In this paper we have studied the nucleoli axial charge and the vector charge in quenched 
lattice QCD. To capture important aspects of the chiral symmetry of QCD, we used domain 
wall fermions to simulate the light quarks. 

We first demonstrated that the lattice renormalization of the isovector vector and axial- 
vector currents satisfy Z v = Z A to a high degree of precision, less than a percent at a^ 1 m 2 
GeV and about two percent at a -1 w 1.3 GeV. This is achieved because in practice the DWF 
method preserves the chiral symmetry of QCD up to small corrections and hence maintains 
the relevant Ward-Takahashi identity. This holds if the underlying (quenched) gauge config- 
uration is sufficiently smooth. Ensembles of such gauge configurations are obtained close to 
the continuum limit. For the single-plaquette Wilson gauge action (3 = 6.0 (a -1 ~ 2 GeV) 
is good enough. For the DBW2 action the lattice spacing may be significantly larger while 
still maintaining good chiral symmetry (a" 1 ps 1.3 GeV). 

Our first calculation of g A with the Wilson gauge action was performed at (3 = 6.0 
on a 16 3 x 32 lattice. The corresponding spatial volume ~ (1.6fm) 3 is similar to those 
used in previous lattice calculations. This volume is rather small in comparison with the 
experimentally measured proton charge radius. On this lattice we found all the relevant 
three-point functions are well behaved and that we can reliably extract the charges. The 
isovector vector current renormalization, Z v , determined from them agrees well with the 
corresponding axial current renormalization, Z A , independently obtained from the axial 
Ward-Takahashi identity. We found that both the axial charge, <?^ ttlce , and its ratio to the 
vector charge, (g A /fiv) lattlce , exhibit a very strong dependence on the quark mass. A simple 
linear extrapolation of (<? A /fiv) lattlcc to zero quark mass yielded a very small value, about 
2/3 of the experimental one. 

The second quenched calculation employed the DBW2 gauge action with (3 = 0.87 set 
for a coarser lattice spacing, a ~ 0.15 fm. This allowed a larger physical volume while 
maintaining good chiral symmetry. To study pure finite volume effects, at fixed lattice 
spacing we calculated on lattices with sizes 16 3 x 32 and 8 3 x 24 (~ (1.2fm) 3 and ~ (2.4fm) 3 , 
respectively) . A significant dependence on the volume is seen in both the axial charge g l ^ tUce 
and the charge ratio {g A /g v ) Uttlcc , with the larger volume giving larger values. In contrast g v 
does not show such dependence. The dependence on the quark mass is also different. In the 
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larger volume the central values remain almost constant, while in the smaller volume they 
decrease noticeably with the quark mass. In the chiral limit the two differ by about 20% 
difference. The behavior of g A on the smaller volume is quite consistent with that observed 
in the earlier calculation with the Wilson gauge action. 

Our estimate of (g A /g v ) Ten at zero quark mass from the larger volume with DBW2 action 
is 1.212 ±0.027(stat) ±0.024(norm). The systematic error is estimated from the two percent 
difference between Z A and Z v which are theoretically equivalent to 0(a 2 ), neglecting even 
smaller effects induced by explicit chiral symmetry breaking in DWF. It underestimates the 
experimental value of 1.2670(30) by less than five percent. This discrepancy is smaller than 
twice the theoretical error. 

Thus dependence on the volume seems to be the largest among the known sources of 
systematic error for the first principles lattice calculation of g A . This suggests that close 
attention be paid to the finite volume effect in other lattice numerical studies of nucleon 
structure, in particular the moments of spin-polarized structure functions which are related 
to the axial charge. 

It should be also noted that although the crucial relation Z v = Z A is satisfied well in 
the second calculation, we detected small 0(a 2 ) differences among different determinations 
of Z v . Such differences were not detectable in the first set of simulations with the Wilson 
gauge action at a 0.10 fm. Numerically this is at most a few percent effect, and does not 
affect the volume dependence. 

As discussed at the end of section III, the present calculation of g A was performed using 
relatively heavy quark masses (390 MeV < < 860 MeV) so that the systematic error 
arising from quenching may be small. However, one may worry that such a calculation 
does not capture the physics of the pion cloud surrounding the nucleon. In spite of this, the 
values of g A n at these unphysically heavy quark masses lie just below the experimental value 
and show little, if any, dependence on the quark mass. It is also interesting to note that 
the product m N g A shows noticeably less volume dependence. Since this is the combination 
that appears in the Goldberger-Treiman relation, it would be interesting to see how its 
counterpart, the induced pseudo-scalar form factor, behaves at small but non-zero nucleon 
momentum. 
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APPENDIX A: CURRENT ALGEBRA AND CVC HYPOTHESIS 



Defining the charge = % J cPxVq (x, t) 4l( , the transformation for the quark fields in 



the isospin SU(2) subgroup of the SU(2) x SU(2) chiral symmetry can be represented by 

[Q a v ,Hx)] = -T>(x), (Al) 
[Q a v 3{x)] = +^(x)T a } (A2) 

where ip = (u,d) T and V£(x) = ip(x) r ) l jT a %p(x). One can easily find the axial current 
A a ^{x) = ^(x)757 M T a -?/'(x) and vector current transform under isospin symmetry as 

[(%X(x)] = te abc V<(x) (A3) 
[Q a v ,Al(x)] = ie^Alix). (A4) 

According to the above SU (2) current algebra, and U+ can be expressed as 

Aj(x) = ^(xh^T^ix) = -[Q$,Al{x)), (A5) 
V?{x) = ^(x) 7m T ± V'(x) = -[Q% Vfc)], (A6) 

where — T\ ± iT 2 . Hence, under the CVC hypothesis, one can find 

(pK\u) = -<p|[GM>> 
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= -{p\Qt A >) + (p\ A lQ» 

= (p\Al\p) - (n\Al\n) 

= 2(p|4|p). (A7) 

The third line follows from Q^\n) = \p) and (p\Q^ — (n\. A similar calculation for the 
vector case yields the relation, (p\V^\n) = 2(p\V£\p). 



22 



[1] K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66, 010001 (2002). 
[2] M. Fukugita, Y. Kuramashi, M. Okawa and A. Ukawa, Phys. Rev. Lett. 75, 2092 (1995) 
[arXiv:hep-lat/9501010]. 

[3] K. F. Liu, S. J. Dong, T. Draper, J. M. Wu and W. Wilcox, Phys. Rev. D 49, 4755 (1994) 
[arXiv:hep-lat/9305025]. 

[4] M. Gockeler, R. Horsley, E. M. Ilgenfritz, H. Perlt, P. Rakow, G. Schierholz and A. Schiller, 

Phys. Rev. D 53, 2317 (1996) [arXiv:hep-lat/9508004]. 
[5] S. Giisken et al. [TXL Collaboration], Phys. Rev. D 59, 114502 (1999). 

[6] D. Dolgov et al. [LHPC collaboration], Phys. Rev. D 66, 034506 (2002) [arXiv:hep- 
lat/0201021]. 

[7] S. Capitani et al, Nucl. Phys. Proc. Suppl. 79, 548 (1999) [arXiv:hep-ph/9905573]. 
[8] R. Horsley [UKQCD Collaboration], Nucl. Phys. Proc. Suppl. 94, 307 (2001) [arXiv:hep- 
lat/0010059]. 

[9] M. L. Goldberger and S. B. Treiman, Phys. Rev. 110, 1178 (1958). 
[10] M. Gell-Mann, R. J. Oakes and B. Renner, Phys. Rev. 175, 2195 (1968). 
[11] D. B. Kaplan, Phys. Lett. B 288, 342 (1992) [arXiv:hep-lat/9206013]. 
[12] Y. Shamir, Nucl. Phys. B 406, 90 (1993) [arXiv:hep-lat/9303005]. 

[13] R. Narayanan and H. Neuberger, Phys. Lett. B 302, 62 (1993) [arXiv:hep-lat/9212019]. 
[14] V. Furman and Y. Shamir, Nucl. Phys. B 439, 54 (1995) [arXiv:hep-lat/9405004]. 
[15] Y. Kikukawa, Nucl. Phys. B 584, 511 (2000) [arXiv:hep-lat/9912056]. 

[16] P. Hernandez, K. Jansen and M. Luscher, Nucl. Phys. B 552, 363 (1999) [arXiv:hep- 
lat/9808010]. 

[17] P. Hernandez, K. Jansen and M. Luscher, arXiv:hep-lat/0007015. 
[18] T. Blum et al, [RBC Collaboration], to appear in PRD, [arXiv:hep-lat/0007038]. 
[19] A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D 63, 114504 (2001) [arXiv:hep- 
lat/0007014]. 

[20] Y. Aoki et al, [RBC Collaboration], [arXiv:hep-lat/02 11023]. 

[21] T. Blum, Nucl. Phys. Proc. Suppl. 73, 167 (1999) [arXiv:hep-lat/9810017]. 

[22] T. Blum et al, Phys. Rev. D 66, 014504 (2002) [arXiv:hep-lat/0102005]. 



23 



T. Blum and S. Sasaki, arXiv:hep-lat/0002019. 

T. Blum, S. Ohta and S. Sasaki, Nucl. Phys. Proc. Suppl. 94, 295 (2001) [arXiv:hep- 
lat/0011011]. 

N. Isgur and G. Karl, Phys. Rev. D 20, 1191 (1979). 

A. Chodos, R. L. Jaffe, K. Johnson and C. B. Thorn, Phys. Rev. D 10, 2599 (1974). 
T. Takaishi, Phys. Rev. D 54, 1050 (1996). 
R. M. Woloshyn and K. F. Liu, Nucl. Phys. B 311, 527 (1989). 
E. Jenkins and A. V. Manohar, Phys. Lett. B 259, 353 (1991). 
G. Martinelli and C. T. Sachrajda, Nucl. Phys. B 316, 355 (1989). 

C. W. Bernard, T. Draper, G. Hockney and A. Soni, in Proc. "Wuppertal 1985 Lattice Gauge 
Theory," 199 (1985). 

M. Kim and S. Kim, Phys. Rev. D 58, 074509 (1998) [arXiv:hep-lat/9608091]. 
S. Sasaki, T. Blum, S. Ohta and K. Orginos [RBC Collaboration], Nucl. Phys. Proc. Suppl. 
106, 302 (2002) [arXiv:hep-lat/0110053]. 
S. Ohta [RBC Collaboration], arXiv:hepdat/0210006. 

L. 1. Wu [RIKEN-BNL-CU collaboration], Nucl. Phys. Proc. Suppl. 83, 224 (2000) [arXiv:hep- 
lat/9909117]. 

R. L. Jaffe, Phys. Lett. B 529, 105 (2002) [arXiv:hep-ph/0108015]. 
T. D. Cohen, Phys. Lett. B 529, 50 (2002) [arXiv:hep-lat/0112014]. 

W. Detmold, W. Melnitchouk and A. W. Thomas, Phys. Rev. D 66, 054501 (2002) [arXiv:hep- 
lat/0206001]. 

S. Sasaki, T. Blum and S. Ohta, Phys. Rev. D 65, 074503 (2002) [arXiv:hepdat/0102010]. 
Note that the Particle Data Group defines g A to be negative because no assumption about 
the structure of the weak interaction is made. In this article, assuming the V — A structure 
of the weak interaction, the axial form factor in Eq. [21 is defined to make g A positive. 
[41] A definition of is slightly different from the one in 3l]]. Only difference is a factor of |. 



24 



TABLE I: Previous lattice calculations with Wilson or O(a) improved Wilson (clover) fermions. In 
general, g A is significantly underestimated. Note that almost all volumes are small, most estimates 
use perturbative renormalization, and partially-unquenching did not increase the value of g A . The 
f symbol denote the continuum extrapolated value. 



type 


group 


fermion 


L 3 x N t 


P 


volume 


statistics 


rn-^L 


g A 


Ref. 




KEK 


V V 1J.DW-L1 


16 3 x 20 


5.7 


(2.2fm) 3 


260 


> 5.9 


985(25") 


[21 




TCentnclcv 




16 3 x 24 


6.0 


(1.5fm) 3 


24 


> 5.8 


1 20(40) 


[31 




DESY 




16 3 x 32 


6.0 


(1.5fm) 3 


1000 


> 4.8 


1 074f90") 


[41 




LHPC-SESAM 

J j 1 1 1 V../ \ — ) A L1V1 




16 3 x 32 


6.0 


(1.5fm) 3 


200 


> 4.8 


1 129(98) 


[61 




QCDSF 


Wilson 


16 3 x 32 


6.0 


(1.5fm) d 


O(500) 




1.14(3)"'' 


[7] 








24 3 x 48 


6.2 


(1.6fm) 3 


O(300) 














32 3 x 48 


6.4 


(1.6fm) 3 


O(100) 










QCDSF-UKQCD 


Clover 


16 3 x 32 


6.0 


(1.5fm) 3 


O(500) 




1.135(34)t 


[8] 








24 3 x 48 


6.2 


(1.6fm) 3 


O(300) 














32 3 x 48 


6.4 


(1.6fm) 3 


O(100) 








Ml(N f = 2) 


LHPC-SESAM 


Wilson 


16 3 x 32 


5.5 


(1.7fm) 3 


100 


> 4.2 


0.914(106) 


[6] 




SESAM 


Wilson 


16 3 x 32 


5.6 


(1.5fm) 3 


200 


> 4.5 


0.907(20) 


[5] 



TABLE II: Simulation parameters for each action and each volume studied in this work. 



Gauge action 


P 


L 3 x N t 


L s 


M 5 


volume 


Wilson 


6.0 


16 3 x 32 


16 


1.8 


(1.6 fm) 3 


DBW2 


0.87 


16 3 x 32 


16 


1.8 


(2.4 fm) 3 






8 3 x 24 


16 


1.8 


(1.2 fm) 3 
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TABLE III: More parameters for each action and each volume. 



Gauge action ((3) 




L 3 x Nt quark mass values statistics (type) m^L 


Wilson (6.0) 




16 3 x 32 0.02 , 0.03, 0.04, 0.05 400 (Wall) >4.3 


DBW2 (0.87) 




16 3 x 32 0.02 , 0.04, 0.06 , 0.08, 0.10 416 (Sequential) >4.8 
8 3 x 24 0.04, 0.06, 0.08, 0.10 400 (Sequential) >3.4 
0.04, 0.06, 0.08, 0.10 400 (Wall) 


TABLE IV: The residual mass m res , hadron masses, inversre lattice spacing (a^ 1 , set by the p 
meson mass) and the renormalization factor of the axial- vector current (Z A ). The p meson mass 
and the nucleon mass are given in the chiral limit in each case. 


Gauge action ((3) M 5 


L s m res nip a p 1 (GeV) Z A Ref. 


Wilson (6.0) 1.8 
DBW2 (0.87) 1.8 


16 1.24 (5) xl0~ 3 0.404 (8) 0.566 (21) 1.922 (40) 0.7555 (3) [18] 
16 5.69 (26) xlO" 4 0.589 (19) 0.780 (27) 1.31 (4) 0.77759 (45) [20] 


TABLE V: Hadron 
L s = 16, from Ref. 


m 
39 


asses computed using Wilson gauge action at /3 = 6.0, 16 3 x 32, M5 = 1.8, 
I- 


m f 




m w m p win 


0.02 
0.03 
0.04 
0.05 




0.2687 (24) 0.4530 (62) 0.645 (12) 
0.3224 (21) 0.4814 (45) 0.716 (5) 
0.3691 (19) 0.5126 (42) 0.754 (6) 
0.4116 (18) 0.5395 (36) 0.805 (5) 


TABLE VI: Results for the nucleon axial charge, Au, Ad and Z v = l/g^ n , Wilson gauge action, 
P = 6.0, 16 3 x 32, M 5 = 1.8, L s = 16, 400 configurations. 


mf (s A )iatt 




(fi^Jrcn (An) rcn (Ad) ren Z v 



0.02 1.216 (106) 0.929 (82) 0.739 (82) -0.189 (40) 0.7637 (23) 

0.03 1.380 (76) 1.056 (59) 0.840 (54) -0.215 (26) 0.7654 (17) 

0.04 1.480 (59) 1.135 (46) 0.903 (40) -0.232 (19) 0.7671 (13) 

0.05 1.542 (47) 1.186 (37) 0.942 (32) -0.243 (15) 0.7689 (11) 
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TABLE VII: Hadron 


masses 


computed using the DBW2 gauj 


$e action. All fits for 


p meson have 


X 2 /Ndf 


< 1.0, and < 1.5 for 


nion and nncleon 






L 3 x N t 




rat 
j 




m p 


m N 






n D4 


0.4255 (38) 


0.679 (16) 


1.071 (13) 


8 3 x 24 




u.uu 


0.5094 (34) 


0.729 (9) 


1.127 (11) 






08 


0.5865 (30) 


0.776 (7) 


1.205 (9) 






0.10 


0.6567 (27) 


0.823 (6) 


1.292 (10) 






0.02 


0.3015 (16) 


0.647 (22) 


0.854 (6) 






0.04 


0.4146 (16) 


0.681 (10) 


0.963 (5) 


16 3 x 32 




0.06 


0.5050 (16) 


0.725 (6) 


1.060 (4) 






0.08 


0.5834 (15) 


0.771 (5) 


1.156 (4) 






0.10 


0.6546 (14) 


0.819 (4) 


1.242 (4) 


TABLE VIII: Results for the nucleoli axial charge, Ait, Ad and Z v = l/g 1 ^, DBW2 


gauge action, 


= 0.87, 8 3 x 24, M 5 


= 1.8, L s = 16, 400 configurations. 






m f 






(# A )ren (Au) rcn 


(AeZ) ren 


Z v 


0.04 


1.303 (146) 


1.059 (120) 0.690 (99) 


-0.369 (99) 


0.8191(65) 


0.06 


1.342 (74) 


1.099 (62) 0.817 (52) 


-0.282 (45) 


0.8242 (35) 


0.08 


1.373 (46) 


1.136 (39) 0.876 (32) 


-0.260 (25) 


0.8317 (24) 


0.10 


1.398 (30) 


1.165 (26) 0.902 (21) 


-0.263 (15) 


0.8403 (18) 



TABLE IX: Results for the nucleon axial charge, Au, Ad and Z y = l/g 1 ^, DBW2 gauge action, 
(3 = 0.87, 16 3 x 32, M 5 = 1.8, L s = 16, 416 configurations. 



m f 




(9 a) ten 


(Af)rcn 


(A<f) rcn 


Zv 


0.02 


1.531 (60) 


1.229 (49) 


0.945 (44) 


-0.284 (27) 


0.8040 (19) 


0.04 


1.523 (24) 


1.230 (20) 


0.946 (17) 


-0.284 (10) 


0.8115 (9) 


0.06 


1.510 (15) 


1.230 (12) 


0.953 (10) 


-0.277 (6) 


0.8184 (6) 


0.08 


1.505 (11) 


1.236 (9) 


0.963 (7) 


-0.273 (4) 


0.8260 (5) 


0.10 


1.503 (8) 


1.246 (7) 


0.975 (6) 


-0.271 (3) 


0.8347 (5) 
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FIG. 1: Z v = l/g l y tt , Wilson gauge action, as a function of the current insertion time-slice. Note 
the very fine scale. A good plateau for each quark mass is observed in the middle range between 
the source and sink. The lines denote central values and statistical errors from constant fits over 
the plateaus. 
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FIG. 2: Quark mass dependence of the vector current renormalization, Z v = l/<jC ltt , Wilson 
gauge action (note scale). Lines denote uncorrelated linear and quadratic fits to the data points. 
Extrapolated values are consistent with the axial-vector current renormalization computed from 

to less than a percent. 
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FIG. 3: The lattice estimates of the spin-dependent densities An (upper curves) and Ad (lower 
curves) for the Wilson gauge action. Decent plateaus are observed for each quark mass. The lines 
denote central values and statistical errors from constant fits over the plateaus. 



30 



2.0 



CD 
O 



CO 



> 

< 



1.5 - 



1.0 - 



0.5 - 



0.0 



0.00 



0.02 



0.04 



0.06 



FIG. 4: The ratio of nucleon charges g A /g v , Wilson gauge action, V ~ (1.6 fm) 3 . A marked decrease 
towards the chiral limit is evident. A linear fit significantly underestimates the experimental value 
(burst). 
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FIG. 5: Z v = l/<?y ltt , DBW2 gauge action, 1/ ~ (2.4 fm) 3 , as a function of the current insertion 
time-slice. Note the very fine scale, same as in Figure ^ We use the three middle points, the 
spread of which is less than 0.5 %. The lines denote the central values and statistical errors from 
constant fits over them. 
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FIG. 6: Same as Figure [2] except for the DBW2 gauge action, large volume 
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FIG. 7: Determination of the vector current renormalization from l/g , |f tt (squares), the electromag- 
netic current (circles), and the d-quark current (diamonds). The axial- vector current renormaliza- 
tion (bursts) 20] is shown for comparison. The different renormalization constants differ because 
of 0(a 2 ) lattice artifacts. 
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FIG. 8: Au and Ad. Same as Figure |3J but for the DBW2 gauge action, sequential source, and 
V = (2.4 fm) 3 . The plateau is shifted towards the wall source because the point sink allows more 
excited state contamination. 
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FIG. 9: Dependence of the unrenormalized nucleon axial charge on quark mass and lattice volume. 
Results from the sequential source method (circles) and the wall source method (bursts) on the 
smaller lattice show good agreement. The sequential method provides somewhat smaller statistical 
errors than the wall method. The larger lattice results (solid circles), obtained with the sequential 
method, exhibit higher values than the smaller lattice ones over the entire range of quark mass 
studied. 
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FIG. 10: The physical ratio of nucleon charges. DBW2 gauge action results on two different 
physical volumes, (2.4 fm) 3 (solid circles) and (1.2 fm) 3 (open circles), reveal the existence of a 
significant finite volume effect. Wilson gauge action results (diamonds), V ~ (1.6 fm) 3 , also appear 
to be affected by finite volume. 
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FIG. 11: Two methods to obtain the physical nucleon axial charge, the ratio of axial- vector to vector 
charge (circles), and the lattice axial- vector charge times the axial- vector current renormalization 



factor Z . in the chiral limit from Ref. 



201 ] . They show slightly different quark mass dependence, 



but exptrapolate to consistent values. Each underestimates the experimental value (burst) by 
rougly five percent. 
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FIG. 12: The lattice volume and quark mass dependences of the product m N g A . Renormalized. 
Closed (large volume) and open (small volume) circles. All agree within one standard deviation. 
No volume dependence is detected. 
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